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Abstract 

Conformal blocks in any number of dimensions depend on two variables z, z. Here 
we study their restrictions to the special "diagonal" kinematics z = z, previously 
found useful as a starting point for the conformal bootstrap analysis. We show that 
conformal blocks on the diagonal satisfy ordinary differential equations, third-order 
for spin zero and fourth-order for the general case. These ODEs determine the blocks 
uniquely and lead to an efficient numerical evaluation algorithm. For equal external 
operator dimensions, we find closed-form solutions in terms of finite sums of 3-F2 
functions. 
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1 Introduction 



Many interesting results about Conformal Field Theory (CFT) in dimensions d ^ 3 have 
been obtained recently using the bootstrap approach [1—15]. Conformal blocks, 1 which sum 
up the contributions of a primary operator and its descendants in the operator product 
expansion, play a fundamental role in this approach, and their detailed knowledge is a 
prerequisite for any advance. Although the theory of these objects was initiated in the 70's 
[16-19], and much developed recently [20-23, 10, 24-27], it is still incomplete. 

To introduce some notation and formulate our problem, consider a correlation function 
of four scalar primary 2 operators fa of dimension Aj. Conformal invariance restricts this 
correlator to have the form 

(M*i)M**)M*s)M**)) = (40 2Al2 (40 2A34 ; rTTTX-TX^TTTX-TT^ . (1-1) 



'12! \^34) 

where x^ = Xi — Xj, Ay = Aj — Aj, and g{u,v) is a function of the conformally invariant 
cross-ratios 

2 2 2 2 

u = ^f, v = ^f. (1.2) 

x 13 x 24 x 13 x 24 



1 In CFT2, one distinguishes the "big" conformal blocks defined summing over the Virasoro descendants, 
and the "small" blocks defined summing over the SL(2, C) descendants only. The CFT^ blocks considered 
here reduce to the "small" blocks for d = 2. 

2 These are called quasi-primaries in CFT2. 
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The four point function (1.1) can be expanded into conformal partial waves correspond- 
ing to the primaries O appearing in the operator product expansions (OPEs) 0i x 2 arid 
03 x 4 . Each conformal partial wave has the same functional form as the four point function 
itself. 3 This gives the following series representation for the function g(u,v): 

g (u, v) = ^2 /120/340 G A/ (u, v) , (1.3) 
o 

Here the GA,e{u,v) are the universal parts of the conformal partial waves; these are the 
conformal blocks mentioned above. They depend on the dimension A and spin I of the 
exchanged primary O, as well as on the external dimension differences A 12 and A 34 . The fijo 
are the OPE coefficients which depend on the CFT under consideration. The representation 
(1.3) has a nonzero radius of convergence [28] and is the starting point for the conformal 
bootstrap studies cited above. 

An equivalent representation of the conformal blocks is expressed in terms of the sym- 
metric functions of variables z, z which are related to u, v via 

u = zz, v = (1 - z)(l - z) . (1.4) 

Then 

G A j(u,v) = F Xl \ 2 (z,z) , (1.5) 

for 

A = Ax + A 2 , £ = Ai-A 2 e {0,1,2,...}. (1.6) 
The normalisation of conformal blocks is a matter of choice, for us 

F\ 1 \ 2 (z 1 z) ~ z A as z — » . (1.7) 



For even d, conformal blocks have closed-form expressions in terms of hypergeometric 
functions [18, 20-22], which proved very useful for conformal bootstrap applications. For 
example, the 4c? blocks are given by [20] 

1 zz 

F Xi\2(z,z)\ = — — r[fc 2 A 1 (^)fc 2 A 2 - 2 (^) -[z+* z)] , (1.8) 

k p (z) = z^ 2 2 F 1 (|(/3 - A 12 ), \(p + A 34 ); P; z) , 

For general d, a numerical approach to computing conformal blocks was developed last year 
in [10], where it was applied in the bootstrap analysis of the 3d Ising model. First, one 
computes the blocks on the "diagonal" z = z, which turns out to be significantly simpler 
than for general z, z. Then, one computes derivatives in the direction normal to the diagonal, 
using recursion relations following from a second-order partial differential equation (PDE) 
that the blocks satisfy. While this method cannot compute the blocks at a finite distance 
from the diagonal, 4 it proved very efficient in computing them in the derivative expansion 

3 See e.g. [23] for an explanation. 

4 See [21, 26] for a proposal on how to do this, using power series expansions around z,z = 0. 
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around z — z — |. The latter information is sufficient for applying the existing conformal 
bootstrap algorithms. 

Motivated by [10], here we analyse in more detail the diagonal limit of conformal blocks. 
The central part of the paper is section 2, where we systematically derive ordinary differential 
equations (ODEs) satisfied on the diagonal. While it is well known that conformal blocks 
as functions of z, z satisfy PDEs (which arise from eigenvalue equations for the Casimir 
operator), it is by no means evident that the diagonal limit satisfies an equation by itself. 
As we will see, this follows from an interplay between the well-known quadratic and the 
rarely-used quartic Casimir PDEs. 

In the rest of the paper we discuss various applications of the ODEs from section 2. 
In section 3 we show that setting one pair of the external dimensions equal, the scalar 
conformal block ODE can be integrated in terms of a 3F2 function, "explaining" the results 
of [10, 29]. For the same case Ai = A2, we also find succinct expressions of higher spin 
blocks as finite sums of 3 F 2 functions. 

In section 4 we propose to use the ODEs from section 2 as a basis for an efficient 
algorithm which can numerically compute conformal blocks and their derivatives around 
any point on the diagonal. This algorithm is a generalisation and an improvement of the 
method used in [10] for the case of equal external scalar dimensions. 

We summarise the main aspects of this paper in section 5. Appendices A and B discuss 
how our ODEs can be used to generate power series expansions of the conformal blocks, 
and a few special cases when the ODEs factorise so that their order can be reduced. 



2 Differential equations on the diagonal 

Conformal blocks in CFT^ may be extended to define symmetric analytic functions of two 
complex variables z and z, everywhere in C 2 except, as follows from (1.8) for d = 4, for 
branch points at z, z = 0,1 and 00. For Euclidean spacetime signature, it is necessary to 
impose z = z* whereas for a Minkowski signature z, z are two independent real variables. 

In this paper we study conformal blocks on the "diagonal" z — z. For the Euclidean 
section z = x + iy,z = x — iy, for x, y real, the limit y — > is in general singular for 
x G [1, 00) 5 . However away from the cut, for arbitrary complex z, the limit z — > z is well 
defined and conformal blocks on the diagonal will be denoted by 

/aiA 2 (z) = F Xl \ 2 (z,z) . (2.1) 

For most conformal bootstrap applications we require real z and < z < 1, which includes 
the crossing symmetric point z — z — \. 

5 This is also apparent from the four-dimensional result (1.8) which can be expressed in the form 
f( z ) g( - z )~^ z ) g ( z ) where f(z), g(z) have branch cuts along the real axis for z > 1. Hence for x > 1, f(z), f(z) 
and g(z),g(z) approach different limits as y — > and the denominator is not cancelled. 
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These functions allow a power series expansion of the form 



f XlX2 (z) =z A J2^nZ n , a = l, (2.2) 



n=0 



where a = 1 follows from the normalisation condition (1.7). 6 Basic CFT properties imply 
z~ f Xl x 2 {z) is analytic in z with singularities at z — 1, oo and so this expansion is convergent 
for \z\ < 1 [28, 26]. 

As functions of z and z, conformal blocks satisfy various PDEs, derived by acting on the 
four point function with the Casimir operators of the conformal group and demanding that 
conformal partial waves be eigenfunctions. For general d, there are quadratic and quartic 
Casimir operators. The corresponding second- [21] and fourth-order [22] PDEs are 

A 2 F Xl \ 2 = c 2 F Xl \ 2 , A 4 F Xl \ 2 = c 4 -FaiA 2 > (2-3) 

where (e = d/2 — 1) 

c 2 = Ai(Ai - 1) + A 2 (A 2 - 1 - 2e) = + 2e) + A(A - 2 - 2e)] , 

C4 = ^ + 2e)(A-l)(A-l-2e), (2.4) 



and the differential operators can be written as 



zz 

A, = D K + D* + 2e 



z — z 



az az 



(2.5) 



A 4 = (j^) £ P- - ^) (^) ' - D M ) , (2.6) 

in terms of the one-dimensional differential operator of the 2 -F\ type: 

d 2 d 

D z = D z {a,b) = (l-z)z 2 — - (a + b+l)z 2 — -abz. (2.7) 

az z az 

Here and below we denote a = — |A 12 , b = |A 34 . 

For £ = the quartic eigenvalue c 4 vanishes. In this case the quartic Casimir equation 
actually reduces to a simple second-order PDE [22] 

(D z -D 2 )F xx = 0. (2.8) 



We will now reduce the PDEs (2.3,2.8) to the diagonal. Conformal blocks are symmetric 
functions of z, z and near the diagonal (2.1) can be extended to a power series expansion of 
the form: 

F XlX2 (z, z) = f XlX2 (t) + \{z - z) 2 g XlX2 (t) + 0((z - z) A ), t = \{z + z) . (2.9) 
6 This is the same normalisation as in [22, 10, 26]. It is related to the normalisation in [20, 21] by 
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The quadratic Casimir equation at z = z gives one relation between / = f\ x \ 2 and g = g\ 1 \ 2 : 



d 2 d 
1(1 -z)z 2 — - (1 + a + b + e)z 2 — - 2ab z - c 2 
dz dz 



f(z) + (l + 2e)(l-z)z 2 g(z) = 0. (2.10) 



A second relation follows from the quartic Casimir equation, which reduces to a third- 
order PDE on the diagonal, so that the 0((z— z) 4 ) terms omitted in (2.9) do not contribute. 
This relation takes the schematic form 

P 3 {d z )f{z) + P 2 {d z )g{z) = 0. (2.11) 

In addition, for I = 0, substituting (2.9) into (2.8) and looking at the 0(z — z) terms gives 
an equation of the form 

Q 2 (d z )f(z) + Q l (d z )g(z) = 0, (2.12) 

The Pi and (J, in the last two equations are certain differential operators of degree i 
with polynomial coefficients, whose precise form is important for what follows but is of no 
particular interest to write them down explicitly. 

With the help of (2.10), we can eliminate g(z) from (2.11) and (2.12). This gives ODEs 
for f(z) by itself, fourth-order for the general case, and third-order for i = 0: 

D^f XlX2 (z) = 0, (2.13a) 
D^f xx (z)=0, (2.13b) 

where in general D^ 1 '^ are differential operators of the form 

D M = {z _ 1) j z «_ + £ (,-ir^(,)^_ + ^ Pr (,)^_, (2.i4) 

r=n—j+l r=0 

for p r (z) polynomials of degree 

I n — r, r > n — j , 
degp r (z) = {. ' (2.15) 
[J, r^n-j. 

The differential operators are symmetric functions of a, b and so, with P = 2ab, S = a + b, 
the operator in (2.13a) is determined by 

p 3 (z) = (AS -2e + 7)z + 4e -2, 

p 2 (z) = [2P + (S - e + 2) (55 - e + 5)}z 2 + 2 [c 2 - P + (3S - e + 3)(2e - l)]z 
+ 4e 2 - 2e - 2c 2 , 

Pi(z) = [4P + (2S + 1)(S - e + 2)](5 - £ + l)z 3 

+ [c 2 (45 - 2e + 1) + 2(25 + 1)(S - e + l)(2e - 1) + P(-AS + lOe - 5)]^ 2 
+ [-6eP + P + 2(25 + l)e(2e - 1) + c 2 (-45 + 6e + l)]z - 2c 2 (2e + 1) , 

p (z) = 2P(S - e)(S -e + l)z 3 + (S - s)[c 2 {2S - 1) + P(Qe - l)}z 2 

+ [c 4 + 2(2e + l)(c 2 (S - 1) + Pe)]z - c 4 + 2c 2 (2e + 1) , (2.16) 
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with c 2 ,c 4 as in (2.4). In (2.13b) 



p 2 (z) = (3S -e + 3)z + 2e, 

Pl (z) = [2P + (2S + 1)(S - e + l)]z 2 + 2(c 2 - P + Se + e)z - 2c 2 , 

p Q (z) = 2P(S - e)z 2 + [P(2e + 1) + c 2 {2S - l)\z + 2c 2 , (2.17) 

where here c 2 = 2A(A — 1 — e). 

With these expressions for the differential operators, the ODEs in (2.13) are a crucial 
result for this paper. They allow for a direct analysis of the diagonal limit of conformal 
blocks. The differential operators satisfy the symmetry relations 

(l-z)^DM {1 - z) -» = D M\^_ ai (2.18) 

for 

' = ^T- < 2 ' 19 > 

and similarly for a <H- b. 

Differential operators of the form (2.14) subject to (2.15) with j ^ n have regular 
singular points at 0, oo, 1 and so are Fuchsian. It is straightforward to see that D^' 
has n(j + 1) — — 1) possible independent parameters. Such differential operators are 
characterised in part by their exponents, which are defined by the roots of the indicial 
equation for series solutions about 0, oo, 1. Thus if Di n '^u(z) = then 

u ~ z so - r , z ->> , u ~ z' Soo - r , z -)• oo, u ~ (z - l) Sl ' r , z-tl, r = 1, . . . , n , (2.20) 

subject to the Fuchs relation 

n 

^(■S ,r + S °o,r + Sl,r) = |n(« - 1) , (2.21) 
r=l 

so that there are a maximum of 3n — 1 independent exponents (although some may be 
degenerate and there may then be logarithmic singularities). For differential operators of 
the form (2.14) with (2.15) and j < n then the exponents for z = 1 are constrained so 
that we may take Si jr = r — 1, r = 1, . . . ,n — j. D^^u^z) = has solutions of the form 
z x n F n ^i(z) with 2n parameters. 

For the differential operators .D^ 4,3 ) and Z)( 3 ' 2 ) of interest here it is easy to determine the 
exponents giving respectively, in terms of A, £, a, b, e, 

oo 1 



A 2a 

2 + 2e-A 2b -2a - 2b (2.22a) 

£ + 1 + 2e a + b — e —a — b — e 

l—£ 1+a+b—e 1—a—b—e 
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and, with I = 0, 

oo 1 

S Ya ° (2 22b) 

2 + 2e-A 2b -2a -2b [ > 

1 a + b — e —a — b — e 

Of course (2.21) is satisfied in both cases. The detailed exponents are useful later in 
discussing possible factorisations of D^' 3 ' ) and Df^ in special cases (see appendix A). 

For general a,b there is no apparent solution of either (2.13a) or more simply (2.13b) 
in terms of known special functions, although results for a = are obtained in section 3. 
However, solutions as a power series expansion around any of the regular singular points 
are generated by Frobenius' method. We are interested in expansions around z = with 
leading term z A . Frobenius' method gives recursion relations by which all the coefficients 
a n in (2.2) can be determined from ao = 1. These recursions, four-term for the general case 
and three-term for I = 0, are given in appendix B. In section 4 we will explain how these 
and related recursions can be used to efficiently evaluate conformal blocks in the derivative 
expansion around any point on the diagonal (and in particular the point z = z = \ relevant 
for the conformal bootstrap applications). 

Mathematically, it is natural to extend the definition of j\ x \. 2 to arbitrary real Ai,A2 
(i.e. I not necessarily a non-negative integer) as solutions of (2.13a). In this case the 
transformations 

A^2 + 2e-A, (2.23a) 
£^ -£-2e, A — >■ A , (2.23b) 

under which both c<i and C4 are invariant, map solutions into solutions. For physical 
dimensions 2e G N. Conformal blocks, with the leading term z Xl+X2 , are related by (2.23a) 
to solutions with leading term /2 2+2£_Al_A2 . On the other hand conformal blocks with the 
same A and spins related by (2.23b) are identical: 

fx t -eX 2 +e(z) =f\ 1 X 2 (z), (2.24) 

since they both satisfy the normalisation condition (2.2). 



3 Closed- form results for a = 

It was found in Ref. [10] that for a = b = the functions fxx(z) and fxx-i(z) can be 
expressed via generalised hypergeo metric functions 3_F 2 . The purpose of this section is to 
give a different derivation of that result, and extend it to higher £. We will also be able to 
generalise this solution to the case when only one of the external dimension differences is 
required to be zero. For definiteness we consider a = and allow b to be arbitrary. 

The essential idea is to consider solutions of the following form: 
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for 



z 2 



Z = z + z' = zz' = ^— Yy (3.2) 

where z' is given by (2.19). Since - — | = ~(i — p-), (3.1) requires that (1 — z) b f(z) is even 
or odd under z -H- z' . 

The utility of the ansatz (3.1) follows from crossing symmetry of the four point function 
under x± <-> x 2 - This implies that conformal blocks satisfy the relation [20]: 

F AlA2 (z', z') = (-lY(l - z)\l - z) b F XlX2 (z, z)\ a ^_ a . (3.3) 

On the diagonal the corresponding relation becomes: 

f XlX2 (z') = e^\l-zf b f XlX2 {z)\ a ^_ a . (3.4) 

The factor e ±l7rA arises since z — )■ z' maps [0, 1) — > (— oo, 0], with ± according which side of 
the branch cut on the negative axis fx^iz') is evaluated on. 

The reason why a = is special is that in this case both sides of (3.4) involve the same 
conformal block and therefore the transformation 

z -> z\ f(z) f(z') = (1 - zf b f{z) (3.5) 

should preserve the defining equations (2.13), as is implied by eqs. (2.18) when a = 0. 
In consequence writing 

f xx (z) = (l-z)- b F X0 (Z), (3.6) 

ensures, since (1— z) b f xx (z) has the same exponents at z = 1, oo and z — > Z maps 1, oo — > oo, 
that we may obtain from (2.13b) a corresponding differential equation for F X0 (Z). This 
equation becomes 

Df 1) F xo (Z) = 0, (3.7) 
involving a simpler D^ 3 ' 1 ^ differential operator: 



= {Z- l)Z 3 ^ + (|(2Z - 1) - e(Z - 1))Z 2 ^ 

- ((6 2 - 1 + e)Z - \c 2 - \e) Z^ + eb 2 Z - \c 2 (3.8) 

= Z(6 Z - e)(S z - b)(6 z + b) - (6 Z - \){5 Z - \)(5 Z + A - 1 - e) , (3.9) 

for 5 Z = Zjg. The form (3.9) leads naturally to a solution in terms of 3F2 functions, as 
mentioned already. Hence the solution of (3.7), consistent with the condition (1.7) at Z = 0, 
is: 
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This result was obtained in [10] for b = (and generalised to nonzero b in [29]) by a different 
method, involving directly manipulating a power series representation of the scalar block. 

Applying the same procedure to (2.13a) leads to an equation involving an operator D^' 2 \ 
This is not manifestly soluble in terms of known special functions. 7 To find the diagonal 
conformal blocks for general I we use instead recursion relations, as discussed in the next 
section. 



3.1 Recursion relations for different I 

Ref. [22] derived differential recursion relations relating F XlX2 (z, z) of different spin (and the 
same a, b). It was noticed in [10] that these relations are particularly useful for relating the 
diagonal limits of conformal blocks, since many derivative terms vanish on the diagonal. 
Here we follow the same idea and set z = z in [22], eq. (4.28), i = 0,1,2. We get the 
following three recursion relations for the f\ 1 \ 2 (z): 

2 ab (c 2 + 2e) \ 

z + 2A 1 (A 1 - 1)(A 2 - e )(A 2 - 1 - £ )) Jx ^ Z} 

" (/aiA 2 -i(z) +*A&i f\ 1+ l\ 2 {z)) 



+ e 



+ (f\i-ix 2 ( z ) +tAp\ 2 -sf\ 1 \ 2 +i(z)) , (3.11a) 

n d , . ab (co + 2e) \ „ 

(1 - .-)- - a - b + (1 + e) 2Ai(Ai _ 1)( l;_ £)( ^_ 1 _ £) ) W*> 

£ + 2e 

" (A 2 /aiA 2 -i(^) -*a/3aj (Ai - 1-£)/ai+1A 2 (^)) 



+ e 



+ — e (( Al + £ ) Ai-iAa^) " *A Px 2 -e (A 2 - 1 - 2e) / Al a 2+1 (z)) , (3.11b) 

(A - 1)(A - 1 - ^ 2A l( A, - iR -l- £ ) 
= (A-l)(/ AlAa _ 1 (z)-/ Al _ 1Aa (z)) 

- (A - 1 - 2e)t A (p Xl f Xl+1X2 (z) - /3 A2 _ £ f Xl X2+1 (z)) , (3.11c) 

where 

(A-2e)(A-l) _ (A + q)(A-a)(A + 6)(A-6) f . 

(A — e)(A — 1 — e) ' 4A 2 (2A - 1)(2A + 1) ' 1 J 

These recurrence relations respect the symmetry (3.4) and also (2.24). They are potentially 
derivable directly from the 4th-order operator Z^ 4,3 ) along the lines of the approach in [22] 
but we do not consider that here. 



7 An alternative D^' 2 ^ equation can be obtained using (3.1) with an extra factor of (1/z — 1/2) on the 
RHS. 
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Eqs. (3.11) may be combined in various ways to remove different terms. Thus 
((A-2-2 £ )(I- 



A-l-2e/ ,\ \\ d 

1 — z a — o 



A-l-e \ v ^ z 2J dz 

+ 2 (A 1 -l)(A 2 -l-e)) /AlA2 (z) 
= ^ (( 2Al " ^ /ax-ia 2 W + (2A 2 - 1 - 2e)(€ + 2e) / Al ^(z)) , (3.13) 
which may be used to determine fw-i starting from £ = 0. 

3.2 Solving the recursion relations 

For conformal blocks for a = and £ = 0, 1, 2, . . . we extend (3.6) to define 

fx\-2n{z)\ a=0 = ^_^ b F X2n {Z) , 
/AA-an-lWU = 773^ (~ - 5) • (3.14) 

The normalisation condition (1.7) requires 

F X2n (Z),F X2n+1 (Z) ~ (-4Z) A ~ n as Z^O. (3.15) 

The recursion relations (3.11) for a = can be applied to obtain corresponding relations 
for F\ n relating even and odd n. From (3.13) for I = 2n we obtain 

2X-2n- 1 -2e . Wl \ 

— — 2n + £ (A - n - 1 - e + 6 Z ) F x 2n (Z) 

IX — In — 1 — s 

= (2A - 4n - 1 - 2e)(n + e) F A2n+1 (Z) + (2A - l)nF A _ 12n _ x (Z) . (3.16) 

Also we may obtain 

(l(2n-l) + 5 z )F X2n (Z) 
= §(2A-l)F A „ 12n _ 1 (Z) 

(2A-2n-l)(A-n-e) (A — 2n — e) 2 - b 2 



(2A - 2n - 1 - e)(2A - 2n - e) 4(2A - 4n+ 1 - 2e) 

which with (3.16) is sufficient to determine F Xn for all n starting from F A o which is given 
by (3.10) 

To solve (3.16) and (3.17) it is convenient to adopt a Mellin transform representation, 
adapting the result for 3 F 2 functions. We therefore assume 



F X2n {Z) = A 



T(2\-2n-e)T{\ + \) 
r(A - n - e) r(A - n + b) T(X-n- b) 
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and 

F A2n+1 (Z) =4 



x 



a-„ r(2A-2n-l-e)r(A + ^) 

r(A - n - e) r(A - n + 6) r(A — n — b) 

27rz J r(s + A — n — 1 — £j r(s + n + |) 



(3.18b) 

The contour of the s-integration is parallel to the imaginary axis between the poles of 
r(— s + A — n) and those generated by the other T functions, moving the contour to the 
right gives a series expansion in powers (— Z) x ~ n+P , p = 0,1,2.... For n = the result 
(3.10) requires 

a xo (s) = l. (3.19) 

The recursion relations (3.16) and (3.17) give 

(A-2n- \-e)(n + s) f3 Xn (s) 
= (A - n - \ - s)(2n + e) a Xn (s) -n(s + n-\) /3 A -in-i(s) , (3.20) 



and 



(2A - 2n - 1 - e) a Xn (s) 
= (s + X-n - 1-e) Px-in-i(s) 

A - n - \ (A - In - ef - b 2 
A-2n + i- e (A - n + 6)(A - n - b) 



+ ( s -X + n) — — j^— -± — — I p Xn .i(s) . (3.21) 



Assuming (3.19), (3.20) and (3.21) then determines a Xn (s) and 0\ n (s) for all n as polyno- 
mials in s of degree n. Clearly from (3.20) j3\o(s) = 1 so that F X i(Z), just as F X0 (Z), is 
given by a single 3F2 function. This generalises eq. (4.11) of [10] to the case of general b. In 
general 

a Xn (\ - n) = (3 Xn {\ - n) = 1 , (3.22) 

as is necessary to ensure (3.15). 

Expressions for a Xn (s), (3 Xn (s) for general n are rather convoluted. Expressions which 
can be extended to arbitrary n can be obtained by considering an expansion in polynomials 
which can be absorbed into the Y functions in (3.18). Then the individual terms in the 
Mellin integral correspond to 3F2 functions with shifted arguments. By trial and error, a 
nice basis is provided by the polynomials 

{s + b) r {s + \) n /{s + \) r , r = 0,...,n. (3.23) 

This ensures T(s + &)—>■ T(s + b + r) and T(s + n + |) _1 — > T(s + r + i)" 1 in the integrand 
of (3.18). Iterating (3.20), (3.21) and expanding in this basis, all the coefficients factorise 
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nicely and the results for a\ n (s), (3\ n (s) become 

{s + \) n (X-e-b-2n) n 
(A - b - ri) n (A - e - 2n + \) n 



a\n\s) 



2 

n i \ _ (s + \) n (X-s-b-2n- l) n 
' Xn{S) ~ {X-b-n) n {X-e-2n-\) n 



s + ~,—\ + e + b + n + l,\ + b — n 



x , { -n,s + b,b + ±,e + n + l | 

1 3 \s + i,-A + e + 6 + n + 2,A + 6-n' j ' "' ' 



Terminating 4^3(1) functions satisfy Whipple identities 



4-T3 



-n, a, 6, c \ _ (e - a) n (f - a) n ^ ( -n,a,d - b,d - c \ 

d,e,f J (e) n (f)n 4 3 \d,a-e + l-n,a- f + l-n /' 

so long as they are balanced, i.e. if 

a + b + c-n + l = d + e + f. (3.27) 

This condition is satisfied in (3.24). The relation (3.26), together with various equivalent 
permutations, then allows (3.24) to be written in a variety of alternative forms, such as 
demonstrating that the results are invariant under b — > —b. Balanced terminating 4^3(1) 
functions also satisfy three term linear contiguous relations, where two of the parameters 
—n, a, b, c, d, e, / are shifted by ±1 while maintaining the balancing condition in (3.27) [30]. 
These are sufficient to show that (3.24) satisfies (3.20) and (3.21). 

With these results in (3.18), the expression (3.10) generalises to a finite sum of 3 F 2 
functions 



Fx2n{Z) = (-4Z) 



x-n (A-n + |) 



X 



E 



(A - n - b) n (A - e - 2n + |) n 
n\ {b + \) r (e + n) r (A - e - b - 2n) n . 
r) (A-n + f) r 



A — n — e, A — n — b, X — n + b + r 
X :; 2 ' 2X-2n- e,X-n + r + \ 



? X2n+1 (Z) = (-4Z)*"" " + t ) 



x 



E 



(A - 7i - fe)„ (A - e - 2n - \) n 
n\ {b + \) r (e + 7i + l) r (A - £ - 6 - 2n - l) n . 



(A-n+|) 



r 



fX-n-s,X-n-b,X-n + b + r \ 
X3F \ 2A-2n-l- £ ,A-n + r + | ' Z ) " (3 ' 28) 

These expressions for higher spin blocks are more concise than those used in Ref. [10]. Using 
eq. (4.9) of [10] the diagonal limit of the spin i ^ 2 conformal block is expressed as a sum 
of i + 1 terms of 3F2 type which is essentially twice the number appearing in (3.28). 
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4 Computing conformal blocks and their derivatives 
efficiently 



We will now return to the case of general a, b, possibly both of them nonzero. In this 
section we will present an efficient algorithm to compute conformal blocks in the derivative 
expansion around any point z = z = to on the diagonal. Such expansions form the basic 
input for the numerical conformal bootstrap algorithms, where to = \ is normally used. Our 
algorithm is an extension of the method first used in [10], with several improvements made 
possible by the results of [26] and of the present paper. 

Let's denote by h mjTl conformal blocks derivatives with respect to the coordinates t, s 
(related by simple rescalings to a, b in [10]): 

F XlX2 (z,z)= V -±-h m , n {t-t G ) m s n , (4.1) 
z — ' mini 

where 

z = t + y/s, z = t-y/s. (4.2) 

The diagonal corresponds to s = and (2.9) corresponds to keeping just n = 0, 1. Conformal 
blocks being symmetric in z z, ensures that the expansion is in integer powers of s. 

The first observation is that the derivatives in the direction orthogonal to the diago- 
nal (s-derivatives) can be recursively determined from the derivatives along the diagonal. 
This recursion follows from the quadratic Casimir equation in (2.3) and has the following 
schematic structure: 

h m ,n= ^ (■■■) h m',n+ [{■ ■ ) h m' ,n-l + {n ~ ■ .)h m > , n ^ 2 ] . (4.3) 

m'^m—1 m'^m+2 

The precise coefficients for a, b = and for t = | were given in [10]; extension to a, b 
nonzero and general < to < 1 is straightforward. By this recursion, moving one unit up in 
n we lose two units in m. So knowing the derivatives h mj0 for m = 0, . . . , m max is sufficient 
to compute h m n for all m + 2n ^ ^ max . In practical applications of numerical bootstrap, it 
is common to use the derivatives in such a triangular table with m max up to O(20) or more. 

Hence, we are reduced to computing the derivatives h m = /i m along the diagonal. We 
next observe that h m satisfy another set of recursion relations as a consequence of the ODEs 
(2.13). These recursions have the schematic form: 

m—1 

m( m -l)(m-2)(m-3)h m = ^ (. . .)h m> (£=1,2,...), (4.4) 

m'=min(0,m— 7) 
m—1 

m(m - l)(m - 2)h m = (. . .)h m > {£ = 0) . (4.5) 

m'=min(0,m— 5) 

The coefficients follow trivially from (2.13) so we do not give them here. Assuming that the 
first few h rn are known (namely for m — 0, 1, 2, 3 for general i and m — 0, 1, 2 for I = 0), 
the rest can be found by these recursions. 
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Thus, it remains to find a method to compute the derivatives h m at low m. This can be 
done as follows. Conformal blocks on the diagonal have an expansion (2.2) around z = 0. As 
we explained in section 2, the expansion coefficients a n are fixed by the ODEs (2.13). The 
closed-form expressions for a n are not available, but they can be found up to an arbitrary 
order via the recursion relations given in appendix B. Crucially, since z = is a regular 
singular point, the single normalisation condition ao = 1 is sufficient to determine all of a n . 
This is unlike the recursions (4.4) for derivatives h n around a regular point < to < 1 where 
derivatives up to the equation order minus one have to be supplied as the initial condition. 

A natural method to compute the h m at low m is then just to evaluate a n up to a 
sufficiently high order N, and to sum up the series (2.2) by differentiating term by term: 

N 

h m = (d z ) m f XlX2 (z)\ z=t0 ^ ^ a„( A + n) ( A + n - 1) . . . ( A + n - m) 4 +n ~ m . (4.6) 

71=0 

We propose to use this method but with a small modification, which greatly improves its 
numerical efficiency. Namely, we will evaluate the conformal blocks and their derivatives 
expanding not in z but in the variable p related to z by 

p= (i + Jr-zr Z = (TVW (47) 

For p the points 0, 1, oo are mapped to 0, 1, —1 and z — >■ z' corresponds to p — > —p. The 
corresponding expansion coefficients will be denoted b n : 

oo 

/ AlA2 (p) = (4p) A ^6 nP ", 6o = l. (4.8) 

ra=0 

The variable p was introduced in [28] as the coordinate corresponding to the OPE frame 
where the operators are inserted symmetrically with respect to the radial quantisation origin. 
Ref. [26] explained that p is a much better in expansion parameter for the conformal blocks 
than z, for two reasons. First, p is smaller that z (e.g. p(|) ~ 0.17) for the ranges of physical 
interest and in applications the series (4.8) converges faster. Second, the coefficients a n grow 
as A n for large A, while b n remain bounded in this limit. So for large A more and more 
terms will have to be retained in the ^-series, while the p-series will not suffer from this 
drawback. 

The ODEs (2.13) in the variable p take the following form: 

2VWp) = 0, V 3 f xx (p) = 0, (4.9) 
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1 d 2 



where 

V* = (P ~ 1) V(P + l) 4 ^ + 2(p - 1) V(p + 1) 3 { (2e + 5)p 2 + 8Sp + 2e - l}^ 

- 2(p - l)p 2 (p + l) 2 {[c 2 - (e + 4)(2e + 3)]p 4 + 4[P - 3S(2e + 3)]p 3 

- 2 [20S 2 + 2e 2 + c 2 + 4P + 3s - 5] p 2 + 4[P + 5(3 - 6e)]p - 2e 2 + c 2 + ej-^ 

- 2p(p + l){(2e + 3)[ca - 2(e + l)]p 6 + [12eP + 6P + 8c 2 S - 8S(e + l)(2e + 3)]p 

+ [4{-4(4e + 3)S 2 - 2e 2 + 8P(S - e) + e + 3} - c 2 (2 £ + 5)] p 4 
- 4 [P(165 - lOe + 5) + 2S (8S 2 + 4e 2 + 2c 2 - 5)] p 3 
+ [-2ec 2 + c 2 + 16P(2S - 2e + 1) - 2 (16S 2 + e - l) (2e - 1)] p 2 

+ [2P(6e - 1) + 8S {-2e 2 + e + c 2 )] p + c 2 + 2c 2 e}^ 

+ (1 - p){ [c 2 (4 £ + 2) - c 4 ]p 6 + 2[-C4 + 2c 2 (25 + l)(2e + 1) + 4P £ (2 £ + l)]p 5 
+ [c4 - 16P(-6eS + S + 2e 2 - 3e) + 2c 2 (l6S 2 + 16Se + 8S + 6e - l)] p 
+ 4 [c 4 + 2c 2 (2S + 1)(45 + 2e - 1) + 4P{85 2 + (6 - 4e)S + e(2e - 3)}] p 
+ [c 4 - 16P{-6eS + 5 + e(2e - 3)} + 2c 2 {16S 2 + 8(2e + 1)S + 6e - 1}] p 
+ 2[-C4 + 2c 2 (25 + l)(2e + 1) + 4P £ (2 £ + l)]p - c 4 + c 2 (4 £ + 2)} , 



7 4 

,3 

2 



and 



% = (P - 1) V(P + l) 3 ^ + 2(P - 1)P 2 (P + lf{{e + 3)p 2 + 6Sp + e}^ 

- 2p(p + l){(c 2 - 2e - 3)p 4 - AS(2e + 3)p 3 

- [2c 2 + 2(8S 2 + e) - 1] p 2 - 8S£ p + c 2 + 4Pp(p - l) 2 ) 

- 2(p - 1) { (c 2 (p + l) 2 + 4P) (p 2 + ASp + 1) + 8P(p + e(p - l) 2 ) } . 



d_ 



They imply recursion relations which determine all b n starting from bo = 1. We do not 
present these recursions here, since they are totally analogous to those for the a n given and 
analysed in appendix B. 

Our method for evaluating h m for low m is thus as follows. First evaluate b n up to a 
sufficiently high order N so that the series 

N 

(d p ) m f Xl x 2 (p)\ P=P0 « J> n (A + n) (A + n - 1) . . . (A + n - m)p^"" m , p = p(t ), (4. 10) 

n=0 

give a good approximation to the RHS for all ^ m ^ m , where m is the maximal 
needed derivative order. The accuracy of this approximation can be controlled via the 
asymptotics of the b n coefficients, which can be understood from the recursion relations 
that they satisfy. For the reasons given above, the needed number of terms N in this series 
will be much smaller than in (4.4). Via the inverse change of variables p — > z, the derivatives 
in z can then be expressed as linear combinations of derivatives in p. 
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5 Summary 



In this paper we have hopefully made several additions to the theory of conformal blocks 
for general dimension d which should be of assistance in implementations of the conformal 
bootstrap: 

1. We derived ODEs satisfied by the blocks on the diagonal z — z. Such equations were 
known before only for spins i = 0, 1 and only for the case of equal external scalar 
dimensions Ai = A2, A3 = A4 [10]. 

2. We proposed an algorithm for an efficient numerical evaluation of conformal blocks 
and their derivatives around any point on the diagonal, generalising the method first 
used in [10] for the case of equal external scalar dimensions. 

3. For the partial case Ai = A2 we derived closed-form representations for the conformal 
blocks on the diagonal in terms of 3F2 functions. Such representations were known 
before only for £ = [10, 29], and for I ^ 1 for the case when the second pair of 
dimensions are also equal: A3 = A4 [10]. For i > 1 our representations are more 
compact and explicit than those in [10]. 

The numerical algorithm from section 4 has an immediate application. Starting from [1] 
and until now, all numerical conformal bootstrap studies have focused on correlators with 
equal external dimensions. As explained in [10] in the context of the 3c? Ising model, it would 
be interesting to perform simultaneous analysis of several four point functions, e.g. (aaaa), 
(aaee), and (eeee), where a and e are the lowest-dimension Z 2 -odd and -even operators. 
Conformal blocks needed for the (aaee) correlator involve unequal external dimensions in 
two out of three OPE channels. They can be computed using our algorithm. 

From a more mathematical perspective the ordinary differential equations obtained here 
may allow more insight into that analytic structure of conformal blocks in general dimensions 
and i > 0. When conformal blocks are expressible in terms of ordinary hypergeometric 
functions, as in even dimensions, their monodromy can be derived from the classic 19th 
century results of Riemann. The extension to n F n _i functions, as appear in some of our 
solutions, was done much later [31] and such techniques may be extendible to the 4th- and 
3rd-order differential operators obtained here. 
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A Factorisation 



Differential equations may in some circumstances be reduced to ones of lower degree. Such 
simplifications are in general a non-trivial mathematical problem [32]. Here we consider 
some cases when the 4th-and 3rd-order differential operators factorise. In general differential 
operators of the form (2.14), subject to (2.15), may factorise so that 

jj(n'+n,f+j) _ £){n',j') jj{n,j) _ (A.l) 

If this is possible the exponents for D^ n ^ must be a subset of those for D^ n ' +n ^' + ^ satisfying 
(2.21). In this case solutions of D^ n +nj = satisfy a lower order differential equation 
if their exponents correspond to those for D^ n '^ . We discuss in this appendix some cases 
when this happens for the operators of interest here. 

First we note that the fourth-order operator contains the third-order one when I = 0, 
| Ai=Aa=A = (% - l)z ^ + (a + b - 1 - e)z + 2e + l) , (A.2) 

where D z 4 ' 3 \ Df^ are determined by (2.16), (2.17) respectively. 

There is also a simplification for conformal blocks with A2 = e — a as a consequence of 

£>0M)| =Df»Df> l \ (A.3) 

z IAi=A,A2=e— a z z i \ j 

(2 1) 

where D z is an operator with exponents 

00 1 



X + e-a 2a (A.4) 

1 — A + £ — a a + b — e —a — b — e 



The form of D z is unimportant but 

Lf.D = ( z - l)z 2 ^- + ((3a + b + 1 - e)z + 2(e - a)) A 

dz 2 dz 

+ 2a(a + b - e)z + (A + e - a) (A - e + a + 1) , (A. 5) 

so that the solution for the diagonal conformal blocks, with A = I + e — a, becomes an 
ordinary hypergeometric function 

fxe-a(z) = z e+2£ - 2a F{£ + e + b-a,£;2e + 2e-2a;z). (A.6) 

For a = b = this result was obtained in [26]. For a = these blocks correspond to the 
"leading twist" operators which saturate the unitarity bound for £ ^ 1, while a > is a 
suitable generalisation. 
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Two other special cases arise for d — 1, 3 and I = when the operator given by (2.14) 
and (2.17) factors as 

Df 2) \ £= _, = - ((z - l)z^- + (a + b - \)z + l) (p,(2a, 26) - 4A(A - §)) , (A.7) 

Df>\ =1 _ = z (l) x {2a, 2b) - 4(A - §)(A - 1)) - [z{l - z) A - ( a + 6 - § )* - 1 ) . (A.8) 

Hence, for d — 1 the solutions with the right asymptotics are of the form z 2X times an 
ordinary hypergeometric function [22], while for d = 3 they can be written as integrals 
thereof: 

= ( i ( !^+i/2 f dt t 1 ~ ^ +6_1/2 * 2A ~ 2i? ( 2A - 1 + 2a, 2A - 1 + 26; 4A - 2; t) . 

(A.9) 

B Recursion relations for a n 

Via Frobenius' method, the ODEs (2.13) imply recursion relations for the coefficients a n in 
(2.2), four-term for general £: 

3 

n(A + £ + n - 1)(A - £ + n - 2e - 1)(2A + n - 2e - 2)a n = ^ 7 i>n a„_;, (B.l) 

i=l 

where 

7l n = 3 ra 4 + (45 + 12A - lOe - 19) n 3 

+ [8e 2 - 2(65 + 15A - 22)e - 4c 2 + 2P + 3(2A - 3) (25 + 3A - 5)]n 2 
+ [12A 3 + 125A 2 - 57A 2 - 365A + 90A + 4(25 + 4A - 5)e 2 + 265 

+ P(4A - 6e - 5) - 2(3A - 4) (45 + 5A - 8)e + c 2 (-45 - 8A + 6e + 13) - 47] n 
+ [A(3A - 4e - 10) + 5(4A - 4 £ - 6) + 6e + 9] (A - 1)(A - 2e - 2) 
+ c 4 + P(A - 2e - 1)(2A - 2e - 3) 

+ c 2 [ - lOe + 5(-4A + Ae + 6) + A(-4A + 6e + 13) - ll] , 

7 2 , n = (A + n + 5 - e - 2){ -3n 3 + [ - 55 - 9A + 5(e + 4)]n 2 

+ [2c 2 - 4P - 3(7^ + 15) + 5(-10A + 8e + 21) + A(10e + 40 - 9 A)] n 
+ c 2 (25 + 2A - 5) + P(-4A + 6e + 7) 

- (A - 2) [5(5A - 8e - 11) + A(3A -he- 14) + lie + 17] }, 

7 3in = (A + 2a + n - 3)(A + 26 + n - 3)(A + n + 5 - £ - 3)(A + n + 5 - £ - 2), 

and three-term for i = 0: 

2 

n(A + n - 1)(2A + n-2e-2)a n = J2 6,n a n -i , (B.2) 
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for 



6, 



a, 



2n 3 + 3(5 + 2A-e-3) 

+ [2P + S(6A - 4e - 9) + 7e + A(5A - 4e - 16) + 13] n 
+ i [2P + (A - 2)(25 + A - 2)] (2A - 2e - 3) , 

- (A + 2a + 7i - 2)(A + 2b + n - 2) (A + S + n - e - 2) . 



Apart from a few cases listed below, these recursion relations determine a n for n ^ 1 
starting from do = 1 and imposing a n = for n < 0. The exceptions occur when the factor 
multiplying a n in the LHS vanishes for some n = uq ^ 1. These are the cases when A 
is smaller by the positive integer, no, than the largest characteristic exponent at z = in 
(2.22). Frobenius' method requires a special treatment in such situations. As we will see 
below, in our case potential ambiguities can be resolved using as an extra physical input 
the fact that conformal blocks should be continuous functions of A. 

There are actually only three exceptional cases consistent with the unitarity bounds and 
assuming d ^ 2. We consider them one by one. 

1. £ > 0, A = £ + 2e, n = 1 in (B.l) 

This corresponds to a spin £ > primary higher spin conserved current saturating the 
unitarity bound. Conformal block of such a primary are defined only for a = b = 0, since 
three point functions with scalars of unequal dimension vanish by imposing conservation. 
Thus we should not worry that eq. (B.l) predicts a\ — > oo for a or b nonzero and A — y £ + 2e. 
However for a = b = conformal blocks should be continuous in this limit. 

Indeed, for n = 1 and a = b = 0, 7^1 factorises: 



The offending factor (A — I — 2e) now appears on both sides of the recursion relation 
defining a±. Cancelling this factor gives an equation for ai which is continuous in the full 
range including the unitarity bound. Once a\ is computed according to this prescription, 
the rest of the coefficients follow from the recursion relation unambiguously. 

We note in passing that the diagonal limit of spin t blocks saturating the unitarity bound 
is known in closed form for arbitrary d [26]; see (A. 6) for a = b = 0. 

2. £ = 0, A = e + |, n = 1 in (B.2) 

Again, these conformal blocks can and should be defined by continuity. Namely, one can 
check that for n — 1 and any a and b the RHS of (B.2) factorises, so that the problematic 
factor (2A — 2e — 1) can be cancelled from both sides of the recursion relation. Once a\ is 
defined this way, the rest of the coefficients follow from the recursion. 

3. £ = 0, A = e,n = 2 in (B.2) 

This case corresponds to a scalar field at the unitarity bound, so it is not particularly 
interesting. This conformal block can occur only in the free scalar theory, and only if the 



7i,i 



i(A - I - 2e)(2A -2s- 1)A(A + £) . 



(B.3) 
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external field dimensions differ by the free scalar dimension e. We have in mind the OPE 
: (j) O : x O D 0. One can check that in this case, i.e. for a,b = ±e/2, the free scalar block 
can be defined by continuity exactly as above. Namely the RHS of (B.2) for n = 2 can be 
factored (plugging in a\ computed in the previous recursion step) and the offending factor 
cancelled. For all the other choices of a, b the scalar block diverges when A — > e. 
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